High-Q impurity photon states bounded by a photonic-band-pseudogap 
in an optically-thick photonic-crystal slab 
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We show that, taking a two-dimensional photonic-crystal slab system as an example, surprisingly 
high quality factors (Q) over 10 are achievable, even in the absence of a rigorous photonic-band- 
gap. We find that the density of in-plane Bloch modes can be controlled by creating additional 
photon feedback from a finite-size photonic-crystal boundary that serves as a low-Q resonator. This 
mechanism enables significant reduction in the coupling strength between the bound state and the 
extended Bloch modes by more than a factor of 40. 

PACS numbers: 42.60.Da, 42.70.Qs, 42.50.Pq 



Significant reduction in the radiation rate of a point- 
like emitter can be achieved by setting up mirrors around 
it,[l[ or by employing photonic crystals (PhCs),@, |1[ 
which is a photonic analogue of atomic crystals for elec- 
tron waves. Q It has long been believed that the existence 
of the photonic-band-gap (PBG) is essential to achieving 
a spatially localized high-Q electromagnetic bound state 
using a PhC cavity. Thus, most efforts so far have focused 
on artificial dielectric structures possessing a PBG.[1-1§| 
Donor- or acceptor-like impurity photon states can be 
formed at the location of a crystal defect. @ Such a lo- 
calized state (with small mode volume, V) has drawn 
much attention in the context of cavity quantum electro- 
dynamics (cQED) experiments 13, U| where the use of 
high Q/V cavities are essential to enhancing light- matter 
interaction. 

Due to fabrication related difficulties, three- 
dimensional (3-D) PhCs have been replaced with a 
lower dimensional counterpart relying on index guiding 
in one or two dimensions. Often, this assumes the from of 
a thin dielectric slab 12j , whose thickness (T) is roughly 
equal to half the effective wavelength (T « A/2n e //) 
in order to optimize the size of the in-plane PBG. One 
representative such a design is shown in Fig. [IJa) and 
(b). Note that a 2-D PhC with the triangular-lattice of 
air-holes supports the PBG only for even guided modes, 
where the symmetry is defined with respect to the mirror 
plane at z — 0. However, this incomplete PBG does not 
preclude the possibility of high-Q defect states, because 
the same mirror symmetry ensures that the defect mode 
with the even symmetry will be completely decoupled 
from all odd guided modes. Lowering the dimensionality 
often creates new symmetries that can be exploited, 
making the condition for localized photon states less 
stringent. 

It is well established that an optically- thick PhC slab 
does not s upp ort any PBG for both even and odd 
symmetries flit Il3j ). As an example, in Fig. QJd), we 
present a photonic band structure (w-fc diagram) for even 



Bloch modes in a thick PhC slab with T = 1.731a. The 
formation of the PBG is hindered by the higher-order 
slab modes lying between the 1st and the 5th bands. 
Note that the energy gap defined by the 1st and the 5th 
bands (Awi_s) wn l remain more or less constant as T 
increases and, in the limit of T — > oo, Ao>i_5 w iH a P~ 
proach the PBG of an ideal 2-D PhC[l2j]. Also note that, 
though all bands, E( Bioc/l )(r), shown in Fig.QJd) are mu- 
tually orthogonal, the bound state in a defect region, 
E( oa „)(r), in general, can couple to any of the higher- 
order Bloch modes. However, this coupling strength, 
Kb | ~ \Jd 3 r AeE {cav) ■ E* (Bloch) \^, shouldn't be 
strong, because Er cav j(z) resembles Ei B i och \{z) of the 1st 
band. In this sense, any optically-thick PhC slab bears 
a pseudo-PBG. Similar pseudo-energy gap in solid state 
physics can create resonance states[15|. Here, the pho- 
tonic counterpart can be easily manipulated with high 
precision by means of the mature modern nanofabrica- 
tion technology. 

In this Letter, we will show how k c b can be controlled 
by setting up simple mirror boundaries around the finite 
size PhC resonator [Fig. [lja)]. The additional boundary 
conditions imposed by the mirrors can alter the density 
of Bloch states in the momentum space (k-space). Thus, 
in relation to the k spectrum of E(„„„) [16j. k Ci b can ex- 
perience significant change. We note that the situation 
is analogous to the well-known cQED example of a point 
dipole source (s=s the hexapole mode) in an optical res- 
onator (r2 the low-Q resonator defined by the boundary 
termination) (l~7| . 

We perform fully-vectorial 3-D numerical simulations 
using the finite-difference time-domain method (FDTD) 



to understand the mutual interaction between E 



(cav) 



and 



E(Bioch)- First, we study the energy decay rate (7) of 
the hexapole modejlH shown in Fig. [IJb). The total 
decay rate (7*0*) can be decomposed into the decay rates 
into the horizontal direction (jhorz) and the out-of-plane 
direction {^ ver t)- Then, 7 is translated into Q factor 
through Q = w/7. Thus, l/Qtot = l/Qhorz + l/Qvert[i3- 
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FIG. 1: (a) A 3-D rendering of the photonic-crystal cavity to be studied in this Letter. The number of photonic-crystal hole 
layers (N layer) and the boundary termination (p) are to be varied, a is the lattice constant of the 2-D triangular-lattice 
photonic-crystal with the background air-holes, 7?, of 0.35a. The refractive index of the slab is assumed to be 3.4. The six 
nearest neighbor holes are reduced (Rm = 0.25a) and slightly pushed away (by a scale factor of 1.1) from the defect center, 
(b) | E | 2 distribution of the hexapole mode in a T — 1.731a slab, (c) Cross-sectional views of the first five lowest Bloch modes 
at M-point. These modes all have even symmetries with respect to the plane at z = 0. (d) The photonic band structure for 
all even Bloch modes. The hexapole mode resonance shown in (b) (u n ~ 0.265) is represented by the horizontal line, in which 
uj n is the normalized frequency defined by Lu n = a/ A. 



Note that we set up detection planes for the Poynting 
energy flux (~Ex H*) away from the mirror boundary, 
so that jhorz accounts for the reflection/transmission at 
the PhC-air discontinuity. 

In Fig. [51 we consider the two boundary terminations 
of p = oo and p = 1.2a. For p — oo, Qtot tends to be 
saturated by Qhorz', Qtot approaches 44,000 ~ 48,000 
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FIG. 2: Qtot, Qvert, and Qhorz of the hexapole mode for two 
different boundary terminations, (a) p — oo and (b) p = 1.2a. 
When the mirror boundary effect is weak (p = oo), Qhorz is 
limited by about 60,000 due to the coupling to the in-plane 
Bloch modes. However, in the case of the abrupt termination 
by air, Qhorz can be made over 2,600,000. 



as N layer > 8[13(. In fact, this particular choice of 
the boundary termination p = oo approximately ensures 
negligibly small reflection off the interface. Alternatively, 
we can simulate the transparent boundary condition by 
overlapping the perfectly matched layer in the FDTD 
with the PhC air- holes [19j, which results in the similar 
saturated Q horz of -60,000 (= Q^'l)- 

Usually, the reflection coefficient at the PhC-air 
boundary is not so large as 0.5~0.6[i3- However, dras- 
tically different behaviors can be seen by forming simple 
mirrors with p = 1.2a [Fig.JSJb)]. (i) All Q values strongly 
modulate with N layer, (ii) Qhorz can be brought up 
to a surprisingly high value of ~ 2.6 x 10 6 at N layer 
= 10. This implies that jhorz can be reduced by more 
than a factor of 40 compared with the case of minimal 
reflection (p = oo). (hi) S ince Qhorz Qvert, now the 
hexapole mode emits more photons into the out-of-plane 
direction. Interestingly, all Qs modulate in the same 
fashion; whenever Qhorz is peaked, so is Q ver t- There- 
fore, the abrupt boundary termination does not simply 
redirect the in-plane guided energy into the out-of-plane 
radiation. Rather, there should be a common physical 
principle for the observed enhancement and suppression. 
Moreover, note that the maximum Qtot of 1.4 x 10 5 in 
the present case can be further increased by employing 
a higher Q ver t cavity mode. Also note that Q ver t of the 
hexapole mode is slightly lowered by the use of finite 
space grids (A = a/20) used in the FDTD simulation [2lj|. 

More detailed analyses are performed as we finely tune 
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FIG. 3: (a-c) Top-down views (|E| 2 ) of the hexapole resonances in a T = 1.731a slab as we increase N layer from 4 to 9 
(in the clockwise direction) for three different boundary terminations of p — 0.8, 1.0, and 1.2a, respectively, (d) The hexapole 
resonances in a thin (T = 0.488a) slab for p = oo are presented in the same manner used in (a-c). Note that this thin slab 
possesses the in-plane photonic-band-gap. (e) Qhorz and (f) Qtot for the hexapole resonances shown in (a-d). 



the mirror boundary conditions. Specifically, we tune p 
in the range of 0.8 — 1.2a as a means to control the phase 
shift (ip) upon reflection. Here, we will provide direct 
graphical evidence for the enhancement /suppression of 
the in-plane Bloch modes. In order to visualize the very 
weak near-fields in the outskirt regions, we adopt the 
saturated color scheme in which 1/1000 of the intensity 
maximum is taken as the upper bound. In Fig. [3(a)- (c), 
we present |E( ca „)(rii, z = 0)| 2 of the hexapole mode for 
three different p = 0.8, 1.0, and 1.2a. We take only one- 
sixth area (cut by the two T-K lines) of the hexapole 
mode profile for each N layer, then combined with plots 
of different N layers to create a one image filling upto 
360°. We also provide a similar plot for a PBG-confined 
hexapole mode in Fig.[3jd) , corresponding to a slab thick- 
ness of 0.488a. For quantitative analyses, we also provide 
graphs for Q tot [Fig. [3(f)] and Qhorz [Fig. [3(g)] as well 
as the resonance frequency shift Au> n /w n [Fig. [3(e)]. 

We find that, at (N layer, p) = (9, 0.8a), 
|E(Bio C ?i)( r ||)| 2 is almost completely suppressed, which is 
comparable with that of the PBG-confined case shown 
in Fig. [3(d). On the other hand, at certain combinations 
of (N layer, p), |E (B;oc/l) (r||)| 2 becomes quite strongly 
pronounced; for example, see (7, 0.8a), (8, 1.0a), and 
(7, 1.2a). For small N layer such as at (N layer < 5, 
p = 1.0a), it seems natural to have more intense photon 
tunneling than the equivalent PBG-confined case. How- 



ever, by slightly tuning p, the strong in-plane loss in the 
previous case can be greatly reduced by more than a fac- 
tor of 10, whose resultant Qtot and Qhorz can be higher 
than those of the PBG-confined case[22j . This result may 
find a practical importance in applications where a device 
miniaturization is required while keeping Q reasonably 
high. 

We find that |E( B ; oc/l )(r||)| 2 variations can explain 
peaks and dips observed in Qhorz or Qtot- We also note 
that the more intensified |E( B ; ocft )(r||)| 2 in the outskirt 
region can contribute to the excessive scattering losses 
into the vertical direction, thereby lowering Q ver t- In 
short, what we have shown here is that k c ,b strongly 
depends on the detailed boundary conditions. k c b ob- 
viously depends on the size of the hexagonal boundary, 
which determines the phase thickness of the 2-D cav- 
ity. It is also understandable that p is a very critical pa- 
rameter controlling the density of in-plane Bloch modes, 
as has been seen in many cQED examples of an atom 
and a cavity. Further evidence of this analogy can be 
found in Fig. [3(e) , which reports the fractional frequency 
shift. Even using sufficiently thick PhC mirrors with iV 
layer > 6, the thick-slab cases show noticeable modula- 
tions in Aw„/w n order of ±10~ 5 , while the PBG-confined 
case does not. These energy- level shifts [13] are signatures 
of the coupling between E( ca „)(r) and E/ Block) ( r ) an d can 
be explained in terms of K Ct B [14j | . 
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FIG. 4: (a) In our coupled-mode model, the hexapole mode is 
assumed to couple only to the Y-M wavevectors. For further 
simplification, the solution for the entire system is projected 
into the x direction. The hexapole symmetry ensures the odd 
mirror symmetry with respect to the x = plane, (b) Qhorz 
enhancement factors in Eq. [4] All plots assume tm = 0.6. 



of the partial mirror shown in Fig. |4|a) and consists of 
the Bloch waves with {/3kjv/}- The CMT states that km 
cannot be arbitrarily determined but rather it should be 
connected by the decay rate into the waveguide port such 
that k m = \l^lhorzer lBM ■ 

When the effect of the photon feedback is weak as 
in the case of p = oo, we can set S+m = in Eq. [T] 
Then, |6| 2 decays exponentially with jhorz = \km\ 2 /2 = 
^o/ Qhorz- However, when there is a feedback mecha- 
nism, jhorz is not a constant but varies with 9m and 4>m 
(reflection phase). For this general situation, we need ad- 
ditional set of equations, which can be derived from the 
CMT[27[. For example, the outgoing power (S-m) and 
the incoming power (S+m) are related each other by b 
through 



8- 



M 



-2i0 



•s. 



+M 



a/27/jo 



(2) 



Though the FDTD provides very accurate first- 
principle means to understand the k c- b modulation, we 
have not been quite convinced as to how the simple 
boundary termination (hence the low-Q) can enable such 
large modulations in Qhorz (and Qtot as well). There- 
fore, it would be instructive to develop a simple model in 
the spirit of the coupled mode theory (CMT) [3 S3 • To 
begin, we would like to note that the in-plane hexapole 
mode profile, |E co „(r||)| 2 , can be approximated in terms 
of the three M-point wavevectors of k M1 = (0,l)|k M |, 
k M2 = (V3/2,l/2)|k A /|, and k M3 = (V3/2,-l/2)|k M | 
with |kftf| = 27r/v / 3a. For example, |E( cat ,)(rn)| 2 « 
| sin(ak M i • — sinfak^ : r ||) — sin(cvk A f 3 • r||)| 2 with 
a correction factor a > l[24j. For our hexapole mode, 
a is ~ 1.1 based on the k-space intensity distribu- 
tion, |E (ca „)(k)| 2 [2l(. This is the reason the outskirt re- 
gion of the hexapole mode resembles the M-point Bloch 
modes [25|]. 

To simplify our CMT model, we assume that the 
hexapole mode couples only to the six M-point wavevec- 
tors, /3x(k M1 , k M 2, k M3 , -k M1 , -k M2 , -k M3 ), where 
j3 is a scale factor ~ 0.7 as evidenced from the photonic 
band structure in Fig.QJd). Now we emphasize that these 
M -point wavevectors ({^kjf }) are closed upon reflections 
at the six hexagonal facets. Thus, we may consider the 
whole set of wavevectors, {/3k&f}, as a channel or a port 
in the CMT formulation [Fig. 2(a)]. Under these assump- 
tions, our hexapole mode can be viewed as side-coupled 
to the port (26J. Then, the time evolution of the hexapole 
mode's energy amplitude (b) can be described by 



db 
dt 



ILJQ 



'fvert 1 "^fhorz 



b + kmS+m, 



(i) 



where the last term describes the total incoming power 
as a result of the feedback by the boundary termination. 
We assume that S+m is measured at the inner boundary 



We define the scattering matrix of the partial mirror such 
that S + m = fMe l ^ M S-m- Then, the modified Qhorz 
(= Qhorz) can b e written using the definition of Q (= ojx 
[total energy stored in the resonator] / [total power loss 
into the port]), 



O' 



feed _ feed 

orz ~~ w £2 



M 



S- 



M 



(3) 



In general, ujf eed (the resonance frequency in the presence 
of the feedback) differs from the original ujq as we have 
seen from the FDTD result in Fig. G3e). However, we 
ujq since the fractional change in u) 



can assume uj 



feed 



is much less than 1%. Note that the partial mirror is 



assumed to be lossless such that r M 



M 



1. 



After solving S-m 
pression for the Qho 

Ihorz/ Ihorzh 



and S+m for 6, we obtain the ex- 
enhancement, QlZt/Qhorl (= 



\feed 
' horz 



= 1 + r M - 2r 2 M cos(2<^aj + 46 M ) 

)(«»*) (1 - r 2 M )(1 + r 2 M + 2r M cos(0 M - 26» M )) ' 

■'horz 



(4) 



Fig. [Kb) shows plots of Eq. @] for several 4>m val- 
ues, while rjw is fixed at 0.6. The model also expects 
drastic modulations in Qhorz depending where we lo- 
cate the resonator with respect to the |E| 2 envelop of 
the low-Q Fabry-Perot type resonator. For example, 
when <pM = 1-Ott, Qhorz enhancement is maximized to 
be about 5 at 8m — tt (=effective half wavelength), 
which can be understood considering that |E| 2 = at 
x = 0. Interestingly, tuning 4>m slightly can greatly im- 
prove the enhancement; 4>m — 0.75-7T results in the max- 
imum QhoTz/ Qhorl verv dose to 30. We also find that 
the 4>m tuning inevitably alters the optimal 9m for Q, 
which agrees with the FDTD results in Fig. [3] the small 
tuning in p changes the optimal N layer for Q. 
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It should be noted that we have assumed only one It- 
port and such a large enhancement holds true only for 
that port. However, if the defect state could excite many 
in-plane Bloch waves with different wavevectors thereby 
creating many k-ports, the overall Qhorz enhancement 
contributed from different k-ports will be averaged out 
to be 1 as Eq. [4] expects. For example, multiple k ports 
might be involved if T were too thick (since this would 
allow more higher-order slab modes with many ks) or if 
the boundary termination would not retain the hexapole 
symmetry (since kj/ cannot be conserved upon reflec- 
tions). 

In summary, we study a quasi-bound photon state 
within a PhC that does not posses a rigorous PBG. We 
show that the density of in-plane Bloch modes can be 
controlled by the termination of the finite-size PhC. The 
coupling strength between the bound state and the ex- 
tended Bloch modes can be reduced by more than a factor 
of 40 in the case of a thick PhC slab. The air-suspended 
slab structure assumed in this study may appear im- 
practical. However, low index cladding material can be 
placed underneath as a supporting structure or a metal 
cladding to the sides of the hexagonal boundary can be 
used a clamp (28|, which could also be used as an elec- 
trode for current injection. By removing the thickness 
constraint of T ~ A/2n e //, many unconventional cav- 
ity designs which were previously discarded because they 
cannot support a PBG can now be reconsidered, such as 
the use of a thick-slab for the design of current-injection 
nanolasers(29j|. 
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